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I. INTRODUCTION 



Observations of binary pulsars [y indicate that the individual neutron stars (NS) in such systems have masses very 
close to the Chandrasekhar limit of 1.4 Mq of white dwarfs. Theoretically, the issue of the maximum mass of a NS 
hings strongly on the equation of state (EOS) and the particle interactions at the high density of the center. Models 
of NS with strong kaon condensation ||^ or even a quark nucleus ('strange star' ||^) can have at most 1.5 - 2 Mq, 
which would leave range for lower mass black holes. However, even modest differential rotation Q may easily increase 
the maximum mass 2 Mq of a nonrotating NS to above 3 Mq for a nascient NS in a transient phase of a supernova. 
Moreover, a mass-quadrupole moment Q is also important B for achieving correspondence with numerical results 




In order to model analytically the exterior field of a NS in the framework of general relativity (GR), one needs an 
exact asymptotically flat solution of the Einstein- Maxwell equations (electrovac spacetimes) possessing at least four 
arbitrary physical parameters which are the mass M, angular momentum J, magnetic dipole n and mass-quadrupole 
moment Q. The simplest solution, besides, can be envisaged as axially symmetric, the magnetic field sharing the 
symmetry of the mass and angular momentum distributions; an additional reflection symmetry with respect to the 
equatorial plane which is expected of the self-gravitating objects (see, e.g., |^ and references therein) must be also 
imposed. 

Our paper aims at presenting an exact solution which does satisfy the above requirements and, what is most 
important, is a mathematically very simple solution admitting a representation exclusively in terms of the rational 
functions of spheroidal coordinates (previous effort in this direction only led either to the solution which had no 
reflection symmetry Q or to solutions [p^ which did not permit the rational function representation, and consequently 
could not be written in a concise form). 



II. FOUR-PARAMETER EXACT SOLUTION 

The reported solution has been constructed with the aid of Sibgatullin's method jl^ according to which the 
complex potentials £ and $ satisfying Ernst's equations are defined, for specified axis data e(z) := £{z,p = 0) 
and f{z) :— ^{z,p = 0), by the integrals 

1 e{^)p{a)d<j 1 f{iMa)da 



If' e(Ofi{a)do If' 

TT J_l Vl - O- TT J_ 



TT J -I Vl - O-^ 

The unknown function /i(cr) is to be found from the singular integral equation 

^1 MW[e(0 + g(r;) + 2/(0/(r7)]rf^ 



((7-r)Vr^^ 
with the normalizing condition 

''^ fi{a)da 



= (2) 



-1 



(3) 



where ^ = z + ipu, r] = z + ipr, p and z being the Weyl-Papapetrou cylindrical coordinates and cr, r G [—1,1]; 
e{r]) := e(ff), f{r]) :— f{ff), and the overbar stands for complex conjugation. 
In what follows, the axis data e(z) and f{z) are chosen in the form 



e(z) 
/(^) 



[z - M - ia){z + ib) + d- 5 - ab 
{z + M - \a){z + ih) + d- 5 - ah' 

|M 

{z + M - \a){z + ib) + d - S - ab' 

p,^ - M^^ 1 
M2-(a-6)2' -"4' 



'^:= .^ ^'^'L d:=][M^-ia-bn (4) 



such that the algebraic equation 



eiz) + e{z) + 2f{z)fiz) = (5) 
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will have a pair of distinct roots of multiplicity two. This is a key point for having the rational form of the final 
expressions for £{p, z) and z) after performing the Riemann-Hilbert procedure of the analytic continuation of 
the functions e{z), f{z) into the complex plane [p, z). The resulting expressions for the potentials £{p, z) and <i>(p, z) 
obtained from Eqs. are of the following polynomial form:^ 

£ = {A - 2MB)/ {A + 2MB), ^ = 2i^C/{A + 2MB), 
A = A[{k^x^ - ~cP ~ ik^xy{a ~ b){x^ ~ 1)] 

-(1 - y^)[(a -b){d~d)~ M^b][{a ~ b){y^ + 1) + 4ikxy], 
B = kx{2k^{x^ - 1) + [bia -b) + 2(5] (1 - y^)} 

+iy{2k^b{x^ - 1) - [fc2(a - 6) - M^b - 2aS]{l ~ y^)}, 
C = 2k'^y{x^ - 1) + [25y - ikx{a - - y^) , (6) 

where we have introduced the generalized spheroidal coordinates 

x^-^{r++r^) and y = ^(r+-r_), 

r± Vp^ + {z± W , k:=VdTS. (7) 

The four arbitrary real parameters entering the solution are the total mass M, total angular momentum per unit 
mass a := J/M, magnetic dipole moment /i and mass-quadrupolc moment 

^ = - 4[M2- t-b)2] + + - (3« + b){a - b)' - V] (8) 

of the source. 

The corresponding complete metric is given by the axisymmetric line elemeni]^ 

d.2 = -m - cd^)^ + k'r' [e'^x' ~ y%-^ + + i)(i - y')dv'] , (9) 

in which the metric coefficients /, 7 and uj are the following rational functions of the coordinates x and y (see, e.g., 
10 for detais of SibgatuUin's method): 

f = E/D, e^^ ^E/16k'{x'-yY, u: = {y^ - 1)L/ E, 

E = {4[fc2(2;2 - 1) + 6(1 - + {a- b)[{a - 6)(rf - S) - M%]il - y^}^ 

-16k^{x^ - 1)(1 - y^){{a - b)[k\x^ - y^) + 25y^] + M^by^}^, 
D = {4{k^x^ - Sy^)^ + 2kMx[2k^{x^ - 1) + (26 + ab - b^){l - y^)] + {a - b)[{a - b) 

x{d~S)- M^b]{y'^ - 1) - A(fy + Ay^{2k^{x^ - l)[kx{a - 6) - Mb] 

-2MM(1 - y2) + [(a _ h)[k^ _ 25) - M%]{2kx + M)(l - y^)Y, 
L = ^k^{x^ -\){{a- b)[k'^{x^ - y^) + 2dy^] + M%y^} 

x{kMx[{2kx + M)^ -a^ + b^ - 2y^ {26 + ab - b^)] - 2y^{'i5d - M%'^)} 

-{4[k\x^ - 1) + 5{\ - y2)]2 + (a _ 5)[(a _ 5)(rf _ <5) _ M%]{1 - y"")^} 

X ((1 - y^){2M{2kx + M)[{a -b){d-5)- b{M^ + 25)] - AM%6 

+ (a - b){A5d - M^b^)} - 8k^Mb{kx + M){x^ - 1)) . (10) 



'^All the formulas of this paper have been checked with the aid of the Mathematica 3.0 computer program [ p^ . 
^Throughout the paper, natural units are used in which the gravitational constant and the velocity of light are equal to unity. 
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Special cases 



Eqs. (6) and (10) admit several well-known classical limits: 

1. In the absence of the magnetic field and vanishing arbitrary quadrupole deformation, i.e /i = 5 = 0, only the 
Ernst potential £ survives which is readily recognizable as that of the Tomimatsu-Sato 5 — 2 solution with 
the mass quadrupole Q = -\{M'^ + 3,P/M). 

2. The stationary pure vacuum limit with a non- vanishing quadrupole parameter & is a particular three-parameter 
specialization of the Kinnersley-Chitre solution . 

3. The magnetostatic limit (a = & = 0) is represented by Bonnor's solution for a massive magnetic dipole. For 
this solution, the quadrupole moment is Q = /i^/M — jM^. 

4. Reduction to the Kerr metric p^ ] with total mass M and total angular momentum per unit mass a :— J/M is 
achieved by setting fi — and then formally choosing = ~ M^. The values of M and a remain independent, 
in particular, < can be imposed since in this special limit the complex continuation 6 — s- 16 is admitted. It 
should be stressed, however, that there exist general arguments p9[ | for the interior of the Kerr metric consisting 
of a perfect fluid according to which a large rotational flattening of the body necessarily implies a large absolute 
value of the quadrupole moment. 

5. It is remarkable that the hyperextreme part of our solution corresponding to pure imaginary values of k belongs 
to the Chen-Guo-Ernst family of hyperextreme spacetimes |^ . This branch might represent exterior flelds of 
relativistic disks. Their importance for astrophysics was shown by Bardeen and Wagoner ||2l[| , cf. also [22| . In 
the absence of the electromagnetic field an exact global solution for an infinitesimally thin disk of dust has been 
constructed by Neugebauer and Meinel cf. |p^ . 

Since neutron stars are known to be 'slowly' rotating astrophysical objects |^5| (even at the Kepler frequency 
ujk = \/ GM/R^ ~ 0.5 ms of a millisecond pulsar would the equator rotate only with a speed of vk ~ c/4), therefore, 
it is the subextreme part of the metric (H) which should be used for their description. At the same time, we still 
need to know the location of singularities in our solution to support the physical interpretation we are attributing to 
it. In Fig. 1 we have plotted in coordinates p and z the typical shapes of the infinite redshift surface which one has 
for the real-valued fc, the dots indicating the position of singularities. The two point singularities on the symmetry 
axis (the poles x = 1, y = ±1) belong to the stationary limit surface, while the ring singularity in each case lies at 
the equatorial plane between the symmetry axis and ergosphere; no singularity outside the infinite red shift surface 
arises, and hence the metric (^) is describing the exterior fields of compact objects such as neutron stars. 

We shall conclude the presentation of our solution by writing out Kinnersley's potential /C |2^] the real part of 
which gives the magnetic component of the electromagnetic four potential: 

/C := + iA'^ = - y^)K/{A + 2MB), 

K = 2k^{x^ - l)[2kx + 3M + iy{a - b)] - (a - b)[2Ma - Mb{l - y^) - A\5y] 

+2{2kx + M)[S{1 - y^) + M{M - iby)] + AM6 , (11) 

where is the electric component of a vector potential associated with the dual electromagnetic field tensor. The 
knowledge of is a necessary basis for the investigation of plasma-dynamical effects around neutron stars. In Fig. 2 
the magnetic lines of force are plotted in cylindrical coordinates for two characteristic cases. 



III. MATCHING TO NUMERICAL MODELS OF NEUTRON STARS 



Our exact axisymmetric solution needs to be matched to interior solutions of neutron stars, in order to be realistic. 
Since the junction conditions on the surface of the NS depend very much on the material details such as equation of 
state, conductivity etc., we will restrict ourselves here only to the identification of asymptotically conserved quantities. 
The identification of mass M and angular momentum J of our solution agrees already with the standard parameters 
of asymptotically flat spacetimes in GR and in the astrophysics pS] of NSs. 
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A. Magnetic field 



Normally, the NS' magnetic field B, predicted already in 1964 by Hoyle, Narlikar, and Wheeler and reaching 
high values below the upper limit of i? < 3 x 10^ Tesla, is ignored in numerical studies of rapidly rotating NS. 
More recently, however, axisymmetric solutions of the Einstein-Maxwell equations have been studied po|] using a 
pseudo-spectral method |^l[| involving Chebyshev-Legendre polynomials in terms of maximal slicing quasi-isotropic 
coordinates. Since the magnetic axis is aligned with the rotation axis and only poloidal fields are permitted, this 
numerical work is particularly suited for a comparison of the electrovac spacetimc outside the NS with our exact 
solution: 

For a star close to a sphere and small polar fields B ^ 10^ Tesla, these numerical results are within an error of 10^'^ 
in agreement with Ferraro's solution [3^ ] 

A^^inpofo^sin^e, r > (12) 
lor 

where is the radius of the star, po is the mass density, and /o is the constant value taken by the electric current 
function |0|. 

The equation above enables us to see how the parameter /i of our solution may depend on the parameters of the 
corresponding interior metric for small values of Q. Indeed, introducing the Boyer-Lindquist-like coordinates R and 
9 via the formulas 

kx^R~M, y^cos9, (13) 
we easily find from (11), in the limit R oo, that 

sin^ 9 „ / 1 \ M sin^ 9 „ / 1 



since R has a representation R = r(l + 0{r~^)) in terms of the isotropic coordinate r From (12) and (14) the 
desired relation of fi to the parameters determining the interior of a NS follows immediately. 



B. Quadrupole moment 

An operational way of defining the quadrupole moment of an axially symmetric body in GR has been developed 
by Ryan |^3|. For NS one finds quite generally ^ that the numerical simulations are rather well accounted by the 
simple quadratic relation 

Q = -c(Af,EOS)^, (15) 

where the constant c — c{M, EOS) depends only on the mass M and the equation of state (EOS) for the interior of 
the NS. For NSs of 1.4 Mq the range of this constant is c = 2 to 7.4. It is intriguing that this simple relation quadratic 
in J holds also for fast rotating NS. 

If we reparametrize the arbitrary quadrupole parameter b of our exact solution rather by the dimcnsionless parameter 
A via b — ±\/ + 2aMA — AP we obtain from (ph for the quadrupole moment (/i = 0): 




/ 1 2 aA2(aA2 - MA ± Va^ + 2aMA ^ 
^~~y^2 2(M2 - a2 + a2A2 _ 2aMA) ^ ' 

Our reparametrization is such that for A = we recover the mass-quadrupole moment Q = —J'^/M of the Kerr 
metric. Thus it is worth pointing out that the mass-quadrupole parameter in our solution is also intimately related 
with the angular momentum dipole and octupole moments, and this means that the deformations of the source are 
mainly due to rotation. In a particular case, for instance, when M = 1.4Mq and a = O.625M0, the values of A 
covering the NS range 2 < c < 7.4 are given by the interval 0.986 < A < 2.068. 

In comparison with (^5|), the constant c for our solution depends not only on the mass M, but also on A > 
(M^ — a'^)/2aM which can be adjusted to different EOS. The additional piece depending on angular momentum per 
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unit mass, however, arises only in higher order of A. Thus the quadrupole moment of our exact solution accounts 
rather well to the simple quadratic law ( |l5| ) of NSs. Moreover, in the 'extreme' limit a — > Af we find 



Q 



( 



1 



A(A ± V2A) 
2(A-2) 



) 



(17) 



whence it can be seen that Q can assume arbitrary values for any given value of M, unlike in the case of the extreme 
Kerr metric for which Q = —M^. In this particular limit, the 'NS interval' corresponds to 0.931 < A < 1.584. 

Futher study is needed to exhibit the astrophysical significance of our solution in more detail. In view of a simple 
analytical form of the new electrovacuum metric and clear physical interpretation of the characteristic parameters it 
possesses, it is anticipated that the solution will prove itself suitable for the use in concrete astrophysical applications 
involving neutron stars. 
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